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Abstract 

We re-examine vacuum polarisation of a scalar field in a quasi-local 
volume including the horizon. We find that Hawking radiation rate 
is derived as a pure decay of vacuum due to scalar field interaction 
with classical gravity exactly in the same way as the origin of vacuum 
polarisation effect in Electrodynamics. 



1 Introduction 

It is often questioned whether Hawking radiation [1] is a real physical effect 
or just a theoretical result due to unreal approximations. This is prompted 
by the fact that for astrophysical solar mass black holes, the Hawking tem- 
perature ~ 10- 8 M Sun /M BH ° K is smaller than the CMB 3° K. The black 
holes formed from stellar collapse satisfy the Chandresekhar limit and are 
bigger in mass than the sun. These will therefore not radiate as the ambient 
temperature is higher than their intrinsic Hawking temperature. Therefore 
experimentally it is unlikely that one will detect this effect, unless there are 
primordial black holes, formed from density fluctuations in the early uni- 
verse. Masses of primordial black holes are small and indeed some might be 
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evaporating away at this time. A search for primordial black holes is yet to 
find results. Surprisingly more hope comes from 'analogue' systems where a 
similar effect can be simulated in the Lab. A very famous example of such 
a system is the sonic 'hole' which has apparently been tested in the Lab 
[2] and more recently the curved 'Graphene sheet' [3]. A horizon might be 
simulated [4] (though there is debate whether it can be implemented in prac- 
tice) using the curved 'sheet of carbon atoms' known as graphene, and one 
expects Hawking radiation to emerge from the system. Given the possibility 
that Hawking radiation can be detected in the Lab, one can re-examine the 
physics of the phenomena and describe it using the existing field theoretic 
paradigms. In the initial derivation of Hawking [1], a geometric approxi- 
mation at the time of horizon formation in a collapsing situation was used 
to show a net flux of particles emergent at asymptotic future. Subsequent 
derivations of particle creation in curved space-times included the Unruh ef- 
fect [6] and vacuum polarisation [7] effects. In this note we will concentrate 
on the derivation of Hawking radiation using a vacuum polarisation effect, 
which is most relevant for describing the effect lets say on a graphene sheet. 

Previous computations of the vacuum polarisation of a scalar field in a 
black hole space-time computed the expectation value of the scalar energy 
momentum tensor and usually the answers were infinite [7]. It was then 
'renormalised' using a particular prescription of subtracting the expectation 
value in Boulware vacuum from say the Unruh vacuum, and in the asymptotic 
limit it yielded the Hawking flux. In this paper, we are interested in studying 
the Hawking flux near the horizon. We define the vacuum polarisation as the 
expectation value of a 'quasi-local action' in the vicinity of the horizon. Local 
quantities like the energy momentum tensor are infinite and one employs a 
point splitting method to extract finite answers [8]. Instead of using this, we 
use the 'action' in a finite region of space and obtain the expectation value of 
this. This is motivated from the Brown and York quasi-local energy concept, 
though we don't exactly use their prescription [9]. 

We show that indeed Hawking radiation flux is recovered as the vacuum 
expectation value of this 'quasi- local action' and thus we can interpret the 
Hawking effect as due to the interaction of scalar vacuaa with classical grav- 
ity which has a horizon. The vacuum relevant for this computation is the 
Unruh vacuum and the finiteness of result (not zero) can be attributed to 
the singularity at the horizon. 

In the next section we describe the vacuum polarisation effect for a scalar 
field in the vicinity of the horizon in details and section three is a conclusion 
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with discussions of open problems in this field. 

2 Vacuum decay in a Quasi-Local Region 

The field theory Fock space vacuum is defined as the state in which there 
are 'no' particles. However, as was observed by Dirac, the vacuum could 
also be interpreted of filled 'negative energy states'. The creation 

of a particle leaves a 'hole' in the Dirac sea, interpreted as the anti-particle. 
Thus it is easy to expect that this 'vacuum' will have a non-trivial behavior. 
As field theory fluctuations happens pair creation and annihilation occurs 
instantaneously in vacuum, and the vacuum has 'virtual particles'. This non- 
triviality of vacuum is manifest as the Casimir effect. Presence of an external 
field can polarise such virtual pairs, and using the available energy, the real 
particles and antiparticles can emerge causing detectable effects. This is 
known as vacuum polarisation. The most discussed 'vacuum polarisation' 
occurs in Electrodynamics, where an external electric field causes the vacuum 
to 'decay' [10]. The exact decay rate for this can be computed using the 
vacuum expectation of the S matrix 

| < 0\S\0 > | 2 = exp(- J d 4 x w(x)) = exp(-W) (1) 

w(x) is the rate of pair creation per unit time per unit volume and W = 
J d 4 x w(x). Obviously if W is non-zero and positive it signifies a 'decay'. 
The S-matrix which encodes the interaction, in case of Electrodynamics, 

is 

S = exp(^ J d A xfA„) (2) 

j 1 * is the operator scalar /fermion matter field current and is the vector 
potential corresponding to a constant classical electric field. 
By expanding the exponential in (1) 

21m r 

W = — < 0| / d^xfA^O > . (3) 

This was computed in [10] and it was crucial that the vacuum expecta- 
tion of the interaction term had a 'imaginary' term for the vacuum to decay. 
We can easily extend the above formula to compute the vacuum polarisa- 
tion of scalar /fermions in presence of an external classical gravitational field. 
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As matter couples to the gravitational field through the energy momentum 
tensor; the interaction is 

J d A x^~g <T TV (4) 

where T^ u is the scalar /fermion energy momentum tensor and g^ the metric 
of space-time, and g its determinant. The rate of vacuum decay would be 

W grav = — <0\j d A x^g- <T T^|0 > . (5) 

What makes this computation slightly difficult is the fact that there is no 
unique fock space vacuum in curved space. In ordinary quantum field theory, 
the vacuum is invariant as observed by a set of inertial observers. In case of 
a general covariant theory an observer is allowed to accelerate, and the QFT 
vacuum does not have a invariant meaning. Thus the expectation values will 
differ according to the particular vacuum in which the observation is being 
done. However, physical results can still be obtained in space times with 
Killing vectors which permit special observers, or when asymptotic properties 
can be used to give a invariant interpretation of the quantum amplitudes. 

We are interested in the computation of such an amplitude in the back- 
ground geometry of a black hole. Hawking radiation was observed as a pro- 
cess of particle creation at the horizon, and it was believed that the origin 
of this is in vacuum polarisation. This (5) amplitude has been calculated 
by several physicists before [11] in the background of a black hole metric, 
but we compute a slightly modified version of (5) motivated from the obser- 
vations of [12]. In this paper using Loop Quantum Gravity (LQG) it was 
shown that the horizon decreases in radius due to the time evolution created 
by a quasi-local Hamiltonian at the horizon. This could be also interpreted 
as the vacuum expectation value of 'gravitational semiclassical fluctuations' 
interacting with the classical geometry of the black hole. The decrease in the 
horizon radius agreed precisely with that expected from the radiation of a 
Hawking quanta. What was crucial in the above derivation was that it was 
vacuum expectation value of: A quasi-local energy at the horizon which was 
approximated by taking limits to the horizon surface from both sides of the 
horizon (inside and outside). 

This suggested that we re-examine the previous derivation of scalar fields 
interacting with the classical gravitational background at the horizon, but 
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where the interaction term is not local but like the Brown and York stress 
energy tensor defined over a finite region [9]. We derive the interaction 
as the boundary term of the scalar action on the surface of a quasi-local 
volume which includes the horizon. We then take the limit where the surface 
coincides with the horizon. We study the interaction at the boundary of a 
hollow cylinder, which has an outer boundary at r# = 2GM + 5 and an inner 
boundary at r A = 2GM — 5. The quasi-local stress-energy of [9] however is 
not suitable in the calculation as it lacks the information of the flux normal to 
the boundaries, and this is precisely what we show contributes to the rate of 
particle creation. We thus compute what we name as the 'quasi- local action' 
instead. We simply take the usual scalar action and compute the value of 
the action over a volume of an annular cylinder which straddles the horizon, 
with one surface at r A outside the horizon and another surface tb inside the 
horizon. The scalar action is, 

s = Jd*x^ <r<^<R0 (6) 

In case there is a boundary defined using a normal n M , the boundary term 
can be found upon simple partial integration 

Vasi = [ * d 3 x ^g g^ 09^0 n v = [ d 3 x^g g^j^- [ d 3 x ^g g^j^n 

Jr A Jtb Jta 

(?) 

where = ifid^cf). Therefore the vacuum decay is simply 

W = 2 Im < 0|S quasi |0 > (8) 

To compute the vacuum polarisation amplitude of this quasi-local action 
(7) one has to find an appropriate vacuum. But the question is of course 
which vacuum as there is no unique vacuum in curved space-time. Currently 
there are three well defined vacuaa in use for the Schwarzschild black hole [7] 

(i) The Boulware Vacuum \B > 

defined by requiring normal modes to be positive frequency with respect to 
the Killing vector with respect to which the exterior region is static. 

(ii) The Hartle-Hawking Vacuum \H > 

defined by taking incoming modes to be positive frequency with respect to 
the Kruskal coordinate V, the canonical affine parameter on the future hori- 
zon, and outgoing modes to be positive frequency with respect to U the affine 
parameter on the past horizon. 
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(iii) The Unruh Vacuum \U > 

denned by taking modes that are incoming from X~ to be positive frequency 
with respect to J^, while those that emanate from the past horizon are taken 
to be positive frequency with respect to U, the canonical affine parameter on 
the past horizon. (See Figure (2)). 

Previous computations of vacuum polarisation using the gravitational en- 
ergy momentum tensor uses a renormalisation prescription or a subtraction 
scheme to compute the Hawking flux, e.g. 

< > ren = Lim™ < U\T^\U >-< B\T^\B > (9) 

is the renormalised Unruh vacuum expectation value of the energy mo- 
mentum tensor at infinity. 

The 'interaction term' or the scalar coupling term is derived from the 
scalar action. We thus concentrate on the action and in our computation we 
simply take the expectation value of the 'quasi- local action' in an appropriate 
vacuum. For the purposes of the calculation it is the Unruh vacuum which is 
most appropriate. In a collapsing situation, particles with positive frequency 
wrt the generators of past horizon emerge from behind the horizon [6]. It is 
these modes which contribute to the inner surface of the cylinder. 

Note we place this cylinder in a collapsing situation, in the region where 
the future horizon has been formed, or the 'steady region' . As the apparent 
horizons and the event horizons coincide in these Cauchy slices, one could 
take a 'quasilocal volume' by taking a spherical region of the apparent horizon 
in one time slice r and generate a volume extending to the next time slice 
r + 5t. 

Thus 

W = 2Im[< U\ I ^g»%n u \U > - < U\ I V=g~g^j,n u \U >] (10) 

Jr B Jta 

We take the vb as a point outside the horizon and ta as a point inside the 
horizon. 

At a point outside the horizon r#, the scalar field operator is a linear 
combination of the modes defined outside the horizon [6] 

To obtain the scalar field in the exterior region, the metric is defined in 
Schwarszchild coordinates 

ds 2 = - ( 1 — jdt 2 +(l —J dr 2 + r 2 dVL (11) 
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Figure 1: The different vacua 



where dVL = d6 2 + sin 2 9d(f) 2 . To solve for the scalar normal modes, we define 

? metric 

2GM\ 



r* = r + 2GM In ( — 1 ) , the metric reduces to 



ds' = - 1 



)[dt 2 -dr* 2 ]+r 2 dn (12) 



The scalar field equations in this background are 

1 



_= «9 M [V^g^dv) = (13) 

d 2 <P-dU+ fl-^— JV = (14) 

(Lg^ is the angular differential operator whose eigenf unctions are the spher- 
ical harmonics Yi m (9,4>)) In the Limit r — > 2GM, the wavefunctions which 
are positive frequency wrt t are well defined as there is a Killing vector d t 
in the exterior region. These are easily identified as oc e - luJt e ±lujr * . Using 
the conventional notation of [8] we identify these as u^im labelled using the 
frequency u and angular momentum quantum numbers /, m: However unlike 
the 'boundary conditions' used in [8], we redefine them for this particular 
physical situation. 
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4>( r B) = J2i a rB U ulm + at B Klm] ( 15 ) 
uilm 

U ^ = JaJ^t 6 '^ Ylm (16) 

This particle was a 'outgoing flux' at the surface of the cylinder r B . 

For the boundary at r^, we take the vacuum corresponding to modes of 
the scalar field which are positive frequency wrt the null generators of the past 
horizon, as these are the modes which emerge from behind the horizon even 
in a collapsing case [6]. As the space-time behind the horizon is accessible 
using the Kruskal coordinates U = e -(*- r *)/ 4GM , V = e (*+ r *)/ 4GM such that 

the metric is defined thus 

\3 



ds 2 = - 32(GM) e- r ' 2m dUdV + r 2 dfi 



(18) 
r 

In these coordinates the future horizon comprises of U — surface. In the 
given figure of the collapsing star the region BDC is inside the horizon. If 
we take the modes as those which are positive frequency wrt d/dll [8, 6]. (U 
reverses sign across the horizon) 

We identify these modes as p^im which are such that duPuiim = —^Puiim- 
Thus straddling the horizon the surface has (this is almost null surface) 

' M = Y] [a rA Puim + a\ A vliJ (19) 



uilm 



According to [6], 

^ = \2 S ^Lm.)]V^ ' M " U ^ + C ^ < 20) 

Using this, we find that < U\S qnas i\U > has an imaginary term which is 
finite at the horizon! Thus putting the (15,16,19,20) in (10) we have in the 
limit r A = r B -> 2GM or S -> 



W = Lim<5_ ) .o2 Im 



J d 3 x^gg rr u ui m d r ul lm - J d 3 x v ^gg rr ^ p u i m d r Pli 

Tb uilm r ' A uilm 

(21) 
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Figure 2: Pair Creation 

W = J cfWnfli- £ ^3I^| 2 ( 22 ) 

where T# = l/87rM is the Hawking temperature. This is precisely the rate 
of particle emission expected and predicted in [1]. Thus the flux naturally 
emerges as in a Bose-Einstein spectrum and is of the form as expected at 
the horizon! A fraction of this will emerge at infinity due to scattering from 
the exterior geometry of the black hole. This computation of the quasilocal 
vacuum polarisation amplitude is a new insight in the method of particle 
creation at the horizon. 

3 Discussions 

In the previous section we showed how the expectation value of a 'quasi- 
local' action at the horizon of a black hole shows that the scalar vacuum 
decays into a flux of particles at a thermal temperature. The interesting 
aspect of this is that instead of a renormalisation prescription, the differ- 
ence of the flux at the inner surface and the outersurface of the timelike 
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cylinder emerges naturally giving the expected Hawking flux. One could 
thus interpret this as the total flux created due to vaccum polarisation at 
the horizon. However, this derivation still does not answer the question of 
how gravity reacts to the escaping flux. To study backreaction one has to 
quantise the matter gravity coupled system and observe the behavior of the 
evolution of the space-time. A semiclassical derivation of time evolving black 
hole space-time exists in the framework of loop quantum gravity [12, 13] and 
this requires a non-unitary evolution. Thus we need further work to really 
describe the phenomena of Hawking radiation. The above calculation can 
be easily extended to include Fermion propagation near the horizon. In ex- 
act analogy with the scalar action, we take the action for the fermion field 
/ d A x^f^g e^7 a (<9 M — w M ) ip where are the tetrads and are the cor- 
responding spin connections. Using the solutions for the Dirac modes which 
are positive wrt the affine generator of the past horizon and those positive 
frequency wrt Schwarzschild time, one should be able to derive the boundary 
contributions / rg d 3 x yj—g ip'y^ip n M which yield the Hawking distribution. 
This will be reported elsewhere [14]. 
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